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Exploring  the  Pathogenic  and  Therapeutic  Implications  of  Aberrant  Splicing  in  Breast  Cancer 

W81XWH-08-1-0402 

PRINCIPAL  INVESTIGATORS:  William  D  Foulkes  MB  PhD,  Jun  Zhu  PhD 

CO-INVESTIGATOR:  Jacek  Majewski,  PhD 

Introduction 

In  this  proposal,  we  aimed  to  a)  systematically  monitor  splicing  variant  profiles  in  breast  cancer 
susceptibility  genes  and  b)  explore  the  role  of  alternative  splicing  in  breast  chemotherapy  using  a  global 
strategy.  In  doing  so,  we  hoped  to  identify  and  validate  candidate  splicing  variants  involved  in 
tumorigenesis  using  polony  digital  exon-profiling  and  functional  assays.  This  annual  report  summarizes  our 
progress  in  the  first  year  of  this  two-year  award. 

Body 

Verbatim  -from  S.O.W 


Objective  1:  To  profile  splicing  patterns  of  100  breast  cancer-related  genes  in  30  normal  and  tumor  breast 
cell  lines 


Objective  2:  To  profile  splicing  patterns  of  100  cancer-related  genes  in  lymphocytes  and  breast  tissues 

Objective  3:  Data  analysis  and  validation  to  determine  the  role  of  aberrant  splicing  in  breast  tumourgenesis 
(months  8-12) 

These  objectives  can  be  summarized  in  terms  of  the  expected  Results/Deliverables  from  Year  1:  1) 
Optimizing  the  method  of  bar-code  PCR-sequencing  for  studying  alternative  splicing  at  candidate  loci.  2) 
Identifying  splicing  variants  associated  with  increased  breast  cancer  susceptibility  derived  from  the 
candidate  loci.  Year  2  focuses  more  on  drug  sensitivity,  and  how  splicing  patterns  may  determine  the 
response  to  certain  chemotherapeutic  agents. 

1)  Optimizing  the  method  of  bar-code  PCR-sequencing  for  studying  alternative  splicing  at  candidate  loci 

As  the  key  subject  of  this  synergistic  IDEA  award  was  the  development  of  a  molecular  barcode,  as  a  prelude 
to  large  scale  analysis  of  splicing  patterns  present  in  the  breast  cancer  genome,  we  undertook  a  the  first 
step  to  evaluate  the  "Molecule  Barcode"  technology  to  monitor  alternative  splicing  in  tumor  samples,  we 
have  performed  a  series  of  proof-of-concept  experiment.  The  results  showed  that  the  effective  sequencing 
length  can  be  virtually  extended  from  36  bp  to  300-400  bp.  While  the  results  are  encouraging,  we  observed 
that  the  coverage  is  biased  towards  the  middle  of  reconstructed  molecules.  In  parallel  with  attempts  to 
improve  our  technology,  and  to  better  monitor  splicing  variants  in  tumor  samples,  we  switched  to  an 
emerging  alternative  strategy,  NimbleGen  Sequence-Capture  array  that  can  selectively  enrich  genes  of 
interest  (3-20  Mbps)  for  deep  sequencing  analysis.  The  captured  fragments  were  then  subjected  to  454 
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pyrosequencing  to  obtain  a  longer  read  length  (~400bp),  which  enables  better  identification  of  alternative 
splicing  and  mutation  detection. 

We  had  worked  with  NimbleGen  to  design  a  custom  capture  array  for  detail  characterizing  a  larger  set 
breast  cancer  related  genes.  The  original  design  was  based  on  385K  array  (Nov,  2008)  and  training  date  was 
set  for  early  2009.  Both  Ting  Ni  (post-doctoral  fellow,  funded  by  this  award,  based  at  Duke  University)  and 
Lili  Li  (PhD  candidate,  funded  from  this  award,  based  at  McGill  University)  took  the  training  courses  using 
supplied  385K  capture  array  by  NimbleGen  during  April  7  -  9,  2009.  Average  288-345  fold  of  enrichment 
using  supplied  genomic  DNA  suggested  the  success  of  training.  We  then  use  the  real  samples  of  the  breast 
cancer  cell  line  BT20  and  drug-resistant  BT20  cell  lines  to  practice  by  ourselves.  Real-time  PCR  results 
suggested  substantial  enrichment  of  desired  regions  (average  enrichment  of  267  fold). 

During  our  interactions  with  NimbleGen,  the  2.1M  capture  array  capture  array  became  commercially 
available  with  an  early  release  in  Feb,  2009.  We  were  offered  by  NimbleGen  to  have  a  free  upgrade  from 
385K  array  to  the  most  advanced  2.1M  array.  Given  that  the  new  custom  array  has  approximately  6-fold 
higher  capacity,  we  decided  to  broaden  our  research  using  the  high-density  array.  The  new  custom  array 
was  design  in  early  March,  which  before  the  actual  training  was  taken  place.  After  the  training,  we 
requested  the  2.1M  custom  array  to  be  delivered.  However,  there  was  a  logistic  mistake  from  Nimblegen 
and  our  order  remained  showing  as  385K  array  in  their  ordering  system.  The  issue  had  been  resolved  in 
June,  2009.  The  array  order  is  currently  in  production  department  and  we  expect  to  have  the  array 
delivered  within  next  2-4  weeks. 

While  waiting  for  the  custom  array  to  be  delivered,  we  have  constructed  two  libraries  with  BT20  cell  line  and 
its  derivative  that  resistant  to  dastanib,  an  important  new  tyrosine  kinase  drug  that  has  considerable  activity 
in  leukemia  and  in  some  solid  tumors,  including  breast  cancer.  We  have  worked  out  an  improved  RNA-Seq 
protocol  to  keep  the  strand  information  during  library  construction.  Therefore,  the  library  will  allow  us  not 
only  to  detect  splicing  variants  but  also  antisense  transcripts  in  cell  transcriptome.  We  have  obtained  initial 
information  that  antisense  transcripts  are  broad  phenomena  of  eukaryotic  transcriptomes  and  may  play 
essential  roles  in  mammalian  gene  regulation.  The  innovative  strategy  we  developed,  in  fact,  can  reliably 
capture  this  hidden  layer  of  mammalian  transcriptomes  and  is  expected  to  provide  novel  insights  into 
altered  gene  expression  in  addition  to  alternative  splicing. 

Lastly,  it  is  worth  pointing  out  the  NimbleGen  capture  array  is  directly  coupled  with  454  pyrosequencing 
platform.  To  prepare  further  study  with  large  number  of  samples,  we  also  tested  the  strategy  to  convert  the 
captured  sequences  to  library  amenable  for  Illumina/Solexa  Sequencing  (which  is  approximately  10  times 
more  cost-effective  than  454  pyrosequencing).  This  led  to  a  optimized  protocol  by  which  the  454  library  can 
be  adapted  for  Solexa  sequencing  through  a  circularization  step. 

In  addition  to  the  analysis  of  the  BT20  lines.  In  addition,  through  collaboration  with  the  Breakthrough  Breast 
Cancer  Institute  in  the  UK,  we  have  collected  RNA  samples  from  34  breast  tumor  lines  that  represent  5 
different  types  of  breast  cancers  classified  by  expression  profiles:  luminal  epithelial-like  subtypes  A  and  B, 
basal  epithelial-like,  ERBB2-amplification  associated  and  normal  breast-like  subtypes.  The  fact  that  the 
response  profiles  of  these  lines  to  Dasatinib  have  been  published,  will  allow  us  to  assess  whether  or  not  the 
genetic  alterations  identified  in  a  induced  resistance  model  can  be  applied  to  intrinsic  resistance.  Thus,  as 
soon  as  our  pipeline  is  well-developed  we  will  be  able  to  test  and  validate  candidate  genes  using  this  well- 
characterised  set  of  breast  cancer  cell  lines. 
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A  third  strategy  we  have  used  is  to  sequence  the  RNA  of  a  breast  cancer  cell  line  directly.  Because  the  pace 
of  technological  innovation  is  substantially  faster  than  the  granting  process,  our  original  proposal  to  analyze 
the  splicing  patterns  of  100  candidate  breast  cancer  genes  (using  the  barcode  and/or  capture  array  strategy) 
has  been  augmented  this  approach  with  state-of-the  art  whole  transcriptome  sequencing.  In  other  words,  in 
this  part  of  the  project,  we  aim  to  analyze  not  the  100  "top"  candidate  breast  cancer  susceptibility  genes, 
but  instead  all  ~24,000  genes.  This  project  has  been  led  by  PI  Foulkes  and  co-l  Majewski.  To  date,  we  have 
used  Illumina/Solexa  RNA  sequencing  to  generate  over  40  million,  50  bp  reads  from  the  transcriptome  of 
the  well  known  breast  cancer  cell  line,  HCC1937.  The  reads  were  mapped  to  the  transcriptome  using  ELAND 
and  CASSAVA  software,  and  visualized  using  the  lllumina  Genome  Studio.  Summary  of  the  mapping  is 
presented  in  Table  1  below. 


Table  1  Summary  of  lllumina  Genome  Analyzer  Sequencing  Reads 


Category 


Total  Number  of  Reads 


Non-Mapping 

Repeat  Masked 

Too  Many  N's  in  Sequence 

Aligned  to  Reference  Genome 

Aligned  to  Splice  Junctions 

Total 


1,515,068 

8,789,665 

1,316 

24,701,485 

5,297,332 

40,304,866 


Over  5  million  reads  are  aligned  to  the  splice  junction  libraries,  and  we  are  using  these  reads  to  quantify  and 
characterize  alternative  splicing  events.  However,  one  of  our  main  interests  is  to  identify  novel  splice 
isoforms  (not  present  in  most  reference  splice  libraries).  To  this  end,  we  are  pursuing  three  approaches:  1) 
Mapping  sequencing  reads  to  exhaustive  splice  junction  libraries  representing  all  possible  splice  junctions 
within  each  gene;  2)  monitoring  excess  reads  that  map  to  normally  intronic  sequences  -  this  is  meant  to 
detect  novel  intron  retention  events;  3)  creating  a  library  and  heuristic  algorithms  for  mapping  reads  to  an 
exhaustive  library  joining  all  exons  across  all  of  the  genes  in  the  human  genome  -  this  last  approach  will 
allow  us  to  detect  all  known  and  novel  alternative  splicing  events,  trans-splicing  (across  genes),  as  well  gene 
fusion  events  that  are  known  to  occur  in  many  tumors.  We  discuss  this  last,  most  inclusive  approach 
below. 

Gene  fusions  are  the  result  of  chromosomal  alterations  involving  two  genes.  These  chimeras  may  have 
severe  phenotypic  effects,  such  as  the  well-studied  BCR-ABL1  fusion  protein  implicated  in  chronic 
myelogenous  leukemia  (Shtivelman  et  al.,  1985)  and  TMPRSS2-ERG  found  in  many  cases  of  prostate  cancer 
(Tomlins  et  al.,  2005).  New  efforts  using  high-throughput  sequencing  have  resulted  in  new  discoveries  of 
gene  fusions.  This  has  prompted  interest  in  determining  whether  these  chromosomal  aberrations  may  be 
specific  to  cancer  and  if  they  are,  may  theoretically  serve  as  an  ideal  diagnostic  and  therapeutic  target 
(Prensner  and  Chinnaiyan,  2009). 

To  detect  gene  fusions  using  our  RNA  sequencing  data,  we  followed  a  similar  approach  described  by  Maher 
et  al.  (2009).  Of  a  total  of  40  million  50  bp  reads  sequenced  from  the  HCC1937  cell  line,  approximately  1.5 
million  reads  were  determined  to  be  non-mapping  after  analyzing  the  data  through  the  lllumina  pipeline 
(version  1.3).  Non-mapping  reads  are  those  which  did  not  have  a  unique  alignment  to  the  reference 
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genome.  We  limited  our  gene  fusion  analysis  to  only  these  reads  as  they  may  potentially  characterize 
breakpoints  of  gene  fusions  not  found  in  the  reference  genome. 

Firstly,  we  downloaded  from  the  UCSC  Genome  Database  (hgl8  assembly)  the  set  of  exon  genomic 
sequences  from  all  mature  mRNA  RefSeq  transcripts.  As  the  mRNA  of  gene  fusions  would  typically  involve 
the  fusion  of  exon  sequences  from  two  different  genes,  we  retained  only  the  boundary  sequences  of  each 
exon  (i.e.  49  bp  sequence  from  the  left  boundary  of  an  exon  and  49  bp  sequence  from  the  right  boundary  of 
an  exon),  in  order  to  identify  reads  that  may  span  such  exon-exon  boundaries.  Blat  (Kent,  2002)  was  used  to 
align  the  non-mapping  reads  against  this  list  of  exon  boundary  sequences.  Using  a  set  of  in-house  scripts 
written  in  Python,  we  filtered  the  results  for  alignments  such  that  a  read  was  partially  aligning  to  an  exon 
boundary.  Of  these  partial  alignments,  we  further  identified  whether  its  remaining  unaligned  sequence 
aligned  to  another  exon  region  either  from  the  Blat  analysis  or  by  simple  string  matching  techniques. 

Figures  1  and  2.  Genome  Studio  screenshots,  one  showing  an  example  of  an  alternatively  spliced  gene 
(FOXM1)  and  one  showing  a  non-expressed  gene,  PTEN,  which  is  known  to  be  deleted  in  this  cell  line. 


Ready. 


Figure  1  Screenshot  from  Genome  Studio  of  PTEN,  which  is  known  to  be  not  expressed  in  the  HCC1937  cell  line.  Here,  only  one 
read  is  mapped  to  this  region  and  is  within  intronic  region  of  PTEN. 
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Figure  2  Example  of  an  alternatively  spliced  gene.  Here,  FOXM1  contains  two  isoforms,  one  of  which  has  two  cassette  exons 
spliced  out.  Based  on  the  reads  summary,  it  can  be  seen  that  the  shorter  isoform  is  being  represented. 


In  addition  to  gene  fusion  detection,  our  approach  will  lead  to  the  detection  of  novel  splicing  events  and 
transcript  isoforms.  It  will  be  possible  to  observe  previously  uncharacterized  exon-exon  junctions  from  a 
particular  gene,  which  may  represent  a  transcript  isoform. 

Finally,  as  part  of  this  project,  we  have  focused  on  splice  variants  of  one  particularly  important  breast  cancer 
susceptibility  gene,  BRCA2.  Previous  studies  have  suggested  that  the  variant  NM_000059.3  c.6853A>G, 
BRCA2:  p. 12285V  (also  known  as  c.7081A>G)  is  unlikely  to  be  a  deleterious  allele  and  in  support  of  this,  we 
observed  that  BRCA2:  p. 12285V  co-occurs  in  trans  with  BRCA2 :  c.5946delT  (more  commonly  known  as 
6174delT).  However,  we  also  noted  that  when  compared  with  wild-type  BRCA2,  BRCA2:  p. 12285V  resulted 
in  increased  exclusion  of  exon  12,  leading  to  the  accumulation  of  a  shortened  form  of  BRCA2  protein.  The 
functional  significance  of  this  shortened  isoform  of  BRCA2  is  uncertain.  Allelic  complementation  in  Brca2 
null  mouse  embryonic  stem  cells  was  performed  by  our  collaborator  Shyam  Sharan  (in  his  letter  of 
collaboration  in  the  grant  submission  he  wrote  “\Ne  are  currently  working  on  several  other  variants, 
including  12285V,  which  although  not  a  high-risk  allele,  is  associated  with  increased  exclusion  of  exon  12. 
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Using  our  system,  we  will  be  able  to  test  splice  variants  that  you  discover".)  He  showed  that  both  p. 12285V 
and  an  allele  completely  lacking  exon  12  (BRCA2Aexl2)  were  functionally  indistinguishable  from  wild-type 
BRCA2.  This  was  based  on  their  sensitivity  to  DNA-damaging  agents,  effect  on  irradiation-induced  Rad51 
foci  formation,  homologous  recombination  and  overall  genomic  integrity.  Taken  together,  our  clinical  and 
functional  studies  strongly  suggested  that  exon  12  is  functionally  redundant  and  therefore  missense  variants 
in  this  exon  are  likely  to  be  neutral.  We  believe  that  such  comprehensive  functional  studies  will  be 
important  adjuncts  to  genetic  studies  of  variants. 

Thus,  although  the  original  the  molecular  barcode  strategy  has  not  yet  been  perfected,  we  are  moving 
forward  on  four  fronts  -  as  stated,  1)  the  barcode  methodology  is  in  development;  2)  we  are  working  with 
state-of-the  art  capture  arrays;  3)  we  are  using  the  very  latest  RNA  sequencing  technology  and  4)  We  are 
conducting  comprehensive  analyses  of  existing  splice  site  variants  in  known  breast  cancer  susceptibility 
genes.  As  stated  above,  in  part  3,  we  have  generated  over  5  million  reads  that  have  been  aligned  to  the 
splice  junction  libraries,  and  we  are  using  these  reads  to  quantify  and  characterize  alternative  splicing 
events.  Moreover,  we  will  add  value  to  this  project  by  attempting  to  identify  fusion  proteins,  as  discussed 
above.  In  part  4,  we  have  determined  that  the  BRCA2  isoform  known  as  BRCA2Aexl2  is  not  associated  with 
a  recognizable  phenotype. 

Key  Research  Accomplishments 

•  Publication  in  press  from  the  Foulkes  group 

•  Development  of  novel  approaches  to  analyzing  RNA  Seq  data 

•  Substantial  enrichment  of  desired  regions  on  the  Nimblegen  array  (average  enrichment  of  267-fold)  has 
been  confirmed  by  real-time  PCR 

Reportable  Outcomes 

•  Manuscript  partly  funded  by  Department  of  Defense:  Lili  Li,  Kajal  Biswas,  Laura  Anne  Habib,  Sergey  G. 
Kuznetsov,  Nancy  Hamel,  Tomas  Kirchhoff,  Nora  Wong,  Susan  Armel,  George  Chong,  Steven  A.  Narod, 
Kathleen  Claes,  Kenneth  Offit,  Mark  E  Robson,  Stacey  Stauffer,  Shyam  K.  Sharan  and  William  D.  Foulkes: 
Functional  redundancy  of  exon  12  of  BRCA2  revealed  by  a  comprehensive  analysis  of  the  BRCA2  variant 
c.6853A>G,  p. 12285V.  Human  Mutation,  in  press,  -  this  ms  is  concerned  with  a  breast  cancer 
susceptibility  gene,  BRCA2.  An  isoform  of  BRCA2  known  as  BRCA2Aexl2  exists.  We  have  explored  the 
functional  significance  of  lacking  exon  12  of  BRCA2. 

•  #1  Seminar  (Dr  Zhu):  "NextGen  Sequencing  analysis  of  eukaryotic  transcriptomes",  NHLBI  (July,  2009) 

•  #2  Seminar  (Dr.  Zhu):  "Single-molecule  barcoding  technology  for  RNA-Seq",  Craig  Venter  Institute 
(April,  2009) 

•  #1  Meeting  Presentation  (Dr.  Zhu):  "Strand-specific  transcriptome  sequencing  for  discovery  of  natural 
antisense  transcripts"  Ting  Ni,  Kang  Tu,  Shen  Song,  Zhong  Wang,  Jun  Zhu,  RNA  2009:  the  Fourteenth 
annual  meeting  of  the  RNA  Society  (May,  2009) 

•  #2  Meeting  Presentation  (Lili  Li)  :  "The  Clinical  Importance  of  Weak  Exons  in  BRCA27'  This  was  a 
question  originating  from  the  comprehensive  analysis  of  12285V  in  Exonl2  -  Annual  report  for  the 
Training  Program  of  Innovative  Technologies  In  Multidisciplinary  Health  Research  University  of 
Manitoba,  Winnipeg,  MB,  Canada  (March  23th,  2009) 

•  #3  Meeting  Presentation  (Lili  Li):  "The  Clinical  Importance  of  a  "Weak"  Exon  in  BRCA2 ? 

Human  Genetics  Graduate  Research  Day,  McGill  University,  Montreal,  QC,  Canada  (May29th,  2009) 
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•  Patent  application:  Zhu  J.,  Ni  T.  and  Gao  Y.  (2008)  Barcode  multiplex  PCR  and  single-molecule 
decomposition.  (US  patent  61/131,279). 

•  #1  Grants  awarded:  Zhu,  Foulkes  and  Majewski  awarded  Susan  G  Komen  for  the  Cure  grant  (2008). 

•  #2  Grants  awarded:  Foulkes  and  Majewski  awarded  Canadian  Breast  Cancer  Research  Alliance  IDEA 
grant  to  discover  novel  fusion  proteins  in  BRCAl-related  breast  cancer  (2009) 

•  #3  Grants  awarded:  Zhu  awarded  a  NSF  grant  entitled  "Genome-wide  exploring  miRNA  Network 
motifs"  NSF-0822033  with  a  recent  "creative  extension"  9/1/2008-8/31/2012. 

•  #1  Grants  applied  for:  Majewski  and  Foulkes  have  applied  for  a  joint  Quebec-China  informatics  grant 

•  #2  Grants  applied  for:  Zhu  has  have  applied  for  a  grant  entitled:  "Functional  elements  in  miRNA 
processing",  N HGRI,  R21 

•  Training:  A  Masters  student  will  join  the  McGill  group  in  September  2009. 


Conclusion 


In  this  proposal,  we  have  set  out  to  evaluate  the  importance  of  splicing  for  breast  cancer  biology.  We  are 
employing  a  number  of  state-of-the-art  technologies  -  in  particular,  we  are  developing  the  very  novel 
"Molecule  Barcode"  technology.  As  this  has  been  a  little  slower  than  we  would  have  liked,  in  addition  we 
have  followed  two  other  avenues  -  a  custom-made  capture  array  to  study  splicing  events  in  breast  cancer 
and  RNA  sequencing,  which  is  a  relatively  unbiased  approach  to  the  problem.  Finally,  we  have  carefully 
evaluated  the  significance  of  one  particular  isoform  of  BRCA2,  and  have  shown  that  it  is  not  associated  with 
an  increased  risk  for  breast  cancer.  This  is  important  for  the  families  whose  members  carry  this  variant.  In 
the  future,  we  will  continue  to  focus  on  developing  our  innovative  technology  development  but  will  also 
pursue  alternative  strategies,  as  described  above.  If  necessary,  we  will  also  continue  to  collaborate  with  Dr. 
Shyam  Sharan  in  attempting  to  classify  splicing  variants  in  the  important  breast  cancer  susceptibility  genes, 
BRCA1  and  BRCA2.  The  knowledge  gained  from  this  project  thus  far  could  directly  benefited  women  with 
breast  cancer  and  those  at  risk  of  breast  cancer.  We  hope  that  our  approach  will  lead  to  significant  insights 
into  the  more  general  question  of  the  importance  of  alternative  splicing  in  breast  cancer  biology. 


References 

Kent,  W.  J.  (2002).  BLAT--the  BLAST-like  alignment  tool.  Genome  Res,  12(4),  656-664. 

Maher,  C.  A.,  Kumar-Sinha,  C.,  Cao,  X.,  Kalyana-Sundaram,  S.,  Flan,  B.,  Jing,  X.,  et  al.  (2009).  Transcriptome 
sequencing  to  detect  gene  fusions  in  cancer.  Nature,  458(7234),  97-101. 

Prensner,  J.  R.,  &  Chinnaiyan,  A.  M.  (2009).  Oncogenic  gene  fusions  in  epithelial  carcinomas.  Curr  Opin 
Genet  Dev,  19(1),  82-91. 

Shtivelman,  E.,  Lifshitz,  B.,  Gale,  R.  P.,  &  Canaani,  E.  (1985).  Fused  transcript  of  abl  and  bcr  genes  in  chronic 
myelogenous  leukaemia.  Nature,  315(6020),  550-554. 

Tomlins,  S.  A.,  Rhodes,  D.  R.,  Perner,  S.,  Dhanasekaran,  S.  M.,  Mehra,  R.,  Sun,  X.,  et  al.  (2005).  Recurrent 
Fusion  of  TMPRSS2  and  ETS  Transcription  Factor  Genes  in  Prostate  Cancer.  Science,  310(5748),  644-648. 


7 


Appendices 


See  attached  ms:  Li  et  al,  Human  Mutation,  in  press,  2009 


8 


Human  Mutation 


Human  Mutation 


Functional  redundancy  of  exon  12  of  BRCA2  revealed  by  a 
comprehensive  analysis  of  the  BRCA2  variant  c.6853A>G, 

P.I2285V 


Journal: 

Human  Mutation 

Manuscript  ID: 

hu  mu-2008-0619. R2 

Wiley  -  Manuscript  type: 

Research  Article 

Date  Submitted  by  the 
Author: 

15-Jul-2009 

Complete  List  of  Authors: 

Li,  Lili;  McGill  University,  Cancer  Genetics 

Biswas,  Kajal;  National  Cancer  Institute,  Mouse  Cancer  Genetics 
Program 

Habib,  Laura;  McGill  University,  Program  in  Cancer  Genetics, 
Departments  of  Oncology  and  Human  Genetics 

Kuznetsov,  Sergey;  National  Cancer  Institute,  Mouse  Cancer 

Genetics  Program 

Hamel,  Nancy;  McGill  University,  Research  Institute  of  the  McGill 
University  Health  Centre 

Kirchhoff,  Tomas;  Memorial  Sloan  Kettering  Cancer  Center, 

Medicine 

Wang,  Nora;  Sir  Mortimer  B.  Davis-Jewish  General  Hospital, 
Department  of  Medical  Genetics  and  Cancer  Prevention  Centre 

Armel,  Susan;  University  of  Toronto,  Familial  Breast  and  Ovarian 
Cancer  Clinic,  Princess  Margaret  Hospital 

Chong,  George;  Sir  Mortimer  B.  Davis-Jewish  General  Hospital, 
Department  of  Pathology 

Narod,  Steven;  University  of  Toronto,  Women's  College  Hospital 
Research  Institute 

Claes,  Kathleen;  Ghent  University  Hospital,  Center  for  Medical 
Genetics 

Offit,  Kenneth;  Memorial  Sloan  Kettering  Cancer  Center,  Medicine 
Robson,  Mark;  Memorial  Sloan  Kettering  Cancer  Center,  Medicine 
Sharan,  Shyam;  National  Cancer  Institute,  Mouse  Cancer  Genetics 
Program 

Foulkes,  William;  Sir  Mortimer  B.  Davis-Jewish  General  Hospital, 
Department  of  Medical  Genetics  and  Cancer  Prevention  Centre; 

McGill  University,  Program  in  Cancer  Genetics,  Departments  of 
Oncology  and  Human  Genetics;  McGill  University,  Research 

John  Wiley  &  Sons,  Inc. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Human  Mutation 


Key  Words: 

BRCA2,  unclassified  variants,  exon  splicing  enhancer,  co¬ 
occurrence,  exon  skipping,  in-frame  deletion,  neutral  variant, 

embryonic  stem  cell 

<$>scholaroNE' 

Manuscript  Central 


John  Wiley  &  Sons,  Inc. 


Human  Mutation 


Page  2  of  34 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Functional  redundancy  of  exon  12  of  BRCA2  revealed  by  a  comprehensive 
analysis  of  the  BRCA2  variant  c.6853A>G,  p.I2285V 

Lili  Li*  ,  Kajal  Biswas*  ,  Laura  Anne  Habib  ,  Sergey  G.  Kuznetsov',  Nancy  Hamel  , 
Tomas  Kirchhoff4,  Nora  Wong5,  Susan  Armel6,  George  Chong7,  Steven  A.  Narod8, 
Kathleen  Claes9,  Kenneth  Oflit4,  Mark  E  Robson4,  Stacey  Stauffer',  Shyam  K.  Sharan2 

William  D.  Foulkes1’3’5 

Program  in  Cancer  Genetics,  Departments  of  Oncology  and  Human  Genetics,  McGill 
University,  Montreal,  QC,  Canada;  2Mouse  Cancer  Genetics  Program,  Center  for  Cancer 
Research,  National  Cancer  Institute  at  Frederick,  Frederick,  MD,  USA;  Research  Institute 
of  the  McGill  University  Health  Centre;  4Clinical  Genetics  Service,  Memorial  Sloan 
Kettering  Cancer  Centre,  New  York,  NY  USA;  department  of  Medical  Genetics  and 
Cancer  Prevention  Centre,  Sir  Mortimer  B.  Davis-Jewish  General  Hospital,  Montreal,  QC, 
Canada;  6Familial  Breast  and  Ovarian  Cancer  Clinic,  Princess  Margaret  Hospital, 
University  of  Toronto,  ON,  Canada;  department  of  Pathology,  Sir  Mortimer  B.  Davis- 
Jewish  General  Hospital,  Montreal,  QC,  Canada;  Women’s  College  Hospital  Research 
Institute,  University  of  Toronto,  ON,  Canada;  9Center  for  Medical  Genetics,  Ghent 
University  Hospital,  Belgium 

*  These  authors  have  contributed  equally 

Corresponding  author: 

William  D  Foulkes 
Room  E-757 

Cancer  Prevention  Centre 
SMBD- Jewish  General  Hospital 

3755  Cote  Ste  Catherine  Road,  Montreal,  H3T  1E2,  QC,  Canada 
Tel.  1  514  340  8222  x3965 
Fax  1  514  340  8712 
william.  foulkes@mcgill.  ca 

Acknowledgments: 

We  would  like  to  thank  Drs  Jun  Zhu  and  Jacek  Majewski  and  Mr  Kevin  Ha  for  helpful 
discussions;  Ms  Megan  Harlan  for  her  help  with  the  New  York  families;  Ms  Sabrina  Notte 
for  help  with  manuscript  preparation  and  the  following  funding  agencies  who  made  this 
work  possible:  Canadian  Breast  Cancer  Research  Alliance  IDEA  and  Department  of 
Defense  grants  (WDF);  National  Cancer  Institute,  National  Institutes  of  Health  of  USA 
(SKS) 


John  Wiley  &  Sons,  Inc. 


1 


Page  3  of  34 


Human  Mutation 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Abstract 

Variants  of  Unknown  Significance  (VUS)  in  BRCA1  and  BRCA2  are  common,  and  present 
significant  challenges  for  genetic  counseling.  We  observed  that  BRCA2:  c.6853A>G 
(p.I2285V)  (BIC  name:  7081A>G)  co-occurs  in  trans  with  the  founder  mutation  c.5946delT 
(p.S1982RfsX22)  (BIC  name:  6174delT),  supporting  the  published  classification  of 
p.I2285V  as  a  neutral  variant.  However,  we  also  noted  that  when  compared  with  wild-type 
BRCA2,  p.I2285V  resulted  in  increased  exclusion  of  exon  12.  Functional  assay  using 
allelic  complementation  in  Brca2  null  mouse  embryonic  stem  cells  revealed  that  p. 12285 V, 
an  allele  with  exon  12  deleted  and  wild-type  BRCA2  were  all  pheno typically 
indistinguishable,  as  measured  by  sensitivity  to  DNA-damaging  agents,  effect  on 
irradiation-induced  Rad5 1  foci  formation,  homologous  recombination  and  overall  genomic 
integrity.  An  allele  frequency  study  showed  the  p. 12285 V  variant  was  identified  in  15/722 
(2.1%)  Ashkenazi  Jewish  cases  and  10/475  (2.1%)  ethnically-matched  controls,  odds  ratio: 
0.99  (95%  confidence  interval:  0.44-2.21),  P  =  0.97.  Thus  the  p.I2285V  variant  is  not 
associated  with  an  increased  risk  for  breast  cancer.  Taken  together,  our  clinical  and 
functional  studies  strongly  suggest  that  exon  12  is  functionally  redundant  and  therefore 
missense  variants  in  this  exon  are  likely  to  be  neutral.  Such  comprehensive  functional 
studies  will  be  important  adjuncts  to  genetic  studies  of  variants. 


Key  words:  BRCA2,  unclassified  variants,  co-occurrence,  exon  splicing  enhancer,  exon 
skipping,  in-frame  deletion,  neutral  variant,  Embryonic  Stem  (ES)  cells 
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Introduction 

Variants  of  unknown  significance  (VUS),  also  known  as  unclassified  variants  (UV), 
account  for  a  small,  but  nevertheless  significant,  proportion  of  all  variants  identified  in 
BRCA1  (OMIM#l  13705)  and  BRCA2  (OMIM#600185).  Interpreting  such  variants  pose 
significant  challenges  for  both  clinicians  and  patients.  The  unclassified  variant 
NM_000059.3:  c.6853A>G  (p. 12285 V)  appears  to  be  restricted  to  the  Ashkenazi  Jewish 
population  and  has  been  reported  79  times  in  the  BIC  database  (http://research.nhgri. 
nih.gov/bic/)  (BIC  name:  7081A>G  )  [Mutation  Nomenclature  used  in  this  paper: 
nucleotide  numbering  reflects  cDNA  numbering,  with  +1  as  the  A  of  the  ATG  initiation 
codon  in  the  reference  sequence  (BRCA2:  NM_000059.3)  (www.hgvs.org/mutnomen).  The 
initiation  codon  is  codon  1 .  Traditional  numbering  based  on  the  BIC  mutation  database  is 
also  provided  (http://research.nh gri.nih.gov/bic/)l.  It  is  the  most  frequently  reported  variant 
in  exon  12,  and  was  identified  97  times  in  68,000  tests  perfonned  by  Myriad  Genetics 
Laboratories.  Several  studies  of  B RCA  1/2  mutations  in  women  with  breast  or  ovarian 
cancer  have  identified  p.I2285V  in  affected  women  and  have  designated  p.I2285V  as  a 
VUS  (Shih  et  ah,  2000;  Pal  et  ah,  2005).  We  were  requested  to  offer  genetic  counseling  to 
an  unaffected  member  of  a  family  who  carries  this  variant.  In  this  context,  it  was  reported 
to  us  (by  a  Canadian  clinical  molecular  genetics  laboratory)  to  be  possibly  associated  with 
aberrant  splicing  of  exon  12  of  BRCA2,  resulting  in  a  truncated  protein.  If  exclusion  of 
exon  12  occurs,  the  result  is  an  in-frame  deletion  of  32  amino  acids  (codon  2281-2312). 

This  isoform  is  known  as  BRCA2AE12. 
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Splice  site  mutations  are  a  common  cause  of  exclusion  of  an  exon  and  many  have  been 
identified  in  BRCA1  and  BRCA2  (Frank  et  al.,  2002).  More  recently,  variants  that  do  not 
alter  the  consensus  splice  site  sequences  but  affect  splicing  have  been  identified  (Liu  et  al., 
2001).  These  mutations  result  in  the  interruption  of  exon  splicing  regulators  and  have  been 
increasingly  recognized  as  important  causes  of  exon  exclusion  or  intron  inclusion  (Cartegni 
et  ah,  2002).  Since  the  BRCA2:  p. 12285 V  was  reported  to  be  associated  with  increased 
exclusion  of  exon  12  in  one  family,  we  were  interested  to  establish  if  this  were  the  case  in 
other  individuals  carrying  this  variant,  and  whether  disruption  of  an  ESE  was  a  likely 
mechanism.  Moreover,  we  wondered  if  exclusion  of  exon  12  was  associated  with  an 
identifiable  phenotype. 

Interestingly,  it  had  been  previously  reported  that  breast  cancer  tissues  show  increased 
amount  of  E12  skipping  transcript  (BRCA2AE12)  (Bieche  and  Lidereau,  1999),  but  no 
specific  function  has  been  attributed  to  this  polypeptide.  Notably,  several  cell  lines  such  as 
MDA-MB-468,  MCF-7,  OVCA420,  DOV13  and  SKOV3  (among  others)  have  increased 
expression  of  this  variant.  In  addition,  the  prostate  cancer  cell  line  DU- 145  expresses 
BRCA2  AE12  at  much  greater  levels  than  the  full-length  isoform  (Rauh-Adehnann  et  al., 
2000).  Despite  these  findings,  co-segregation  and  family  history  analyses  have  shown  that 
p.I2285V  is  unlikely  to  be  a  highly  penetrant  allele  (Easton  et  al.,  2007).  This  view  is 
supported  by  the  recently  published  work  of  Tavtigian  and  colleagues,  which  predicts  that 
p.I2285V  is  a  neutral  variant,  because  the  two  independent  measures  used  (missense 
mutations  outside  of  the  DNA  binding  domain  and  family  history  likelihood  ratio)  place 
this  variant  in  a  group  that  has  a  prior  probability  of  <0.05  (Tavtigian  et  al.,  2008).  The 
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potential  effect  of  a  variant  on  splicing,  however,  is  not  measured  in  these  analyses.  Taking 
all  the  observations  together,  it  seems  possible  that  this  variant  affects  splicing  and  thus  has 
a  subtle  effect  on  risk,  perhaps  akin  to  the  risk  associated  with  so-called  low  penetrance 
breast  cancer  genes  identified  recently  (Pharoah  et  ah,  2008).  Here,  we  describe 
investigation  of  the  relationship  between  this  variant  and  exclusion  of  exon  12;  b)  the 
functional  consequences  of  exclusion  of  exon  12  of  BRCA2  and  c)  the  potential  effect  of 
this  variant  on  breast  cancer  risk. 

Materials  and  Methods 
BRCA2:  p.I2285V pedigrees 

We  collected  information  from  four  families  where  at  least  one  individual  carried  the 
BRCA2:  p.I2285V  variant.  These  pedigrees  are  shown  in  supplementary  figures  la-d. 
Lymphocyte  RNA  and  DNA  was  extracted  from  one  carrier  of  the  variant  from  each  family 
and  were  used  in  the  experiments  described  below.  All  participants  signed  an  informed 
consent  form.  The  details  of  each  family  are  described  in  the  figure  legend  accompanying 
supplementary  figure  1 . 

Testing  for  inclusion  or  exclusion  of  exon  12  in  humans 

We  confirmed  the  co-occurrence  of  c.6853A>G  /  p.I2285V  (BIC  name:  7081A>G )  with 
another  deleterious  BRCA2  mutation  c.5946delT  /  p.S982RfsX20  (BIC  name:  6174deIT), 
using  a  blood  sample  from  a  patient  carrying  both  alleles  (see  supplementary  figure  2  and 
accompanying  legend).  We  used  RNA  extracted  from  blood  samples  from  this  proband 
and  three  others  (supplementary  figure  1)  using  PAXgene  (Qiagen)  kit  and  cDNA  was 
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synthesized  using  Sensiscript  (Qiagen).  We  used  the  following  PCR  primers  to  semi- 
quantitatively  compare  the  degree  of  exon  12  skipping  between  the  cells  lines  the  wild-type 
line  (referred  to  as  WLT)  and  the  double  variant  line  (referred  to  as  DBV):  1)  for  exon  12 
inclusive  transcript:  Forward  (located  in  exon  12)  :  5’  -  CCCTTATCTTAGTGGG 
AGAACCCTCA-3’,  Reverse  (located  in  exon  13):  5'- 

AGTTGT  GCGA  AAGGGT  AC  AC  AGGT-3 ' ;  2)  for  constitutive  exon  11:  Forward  (located 
in  exon  1 1):  5’  -TT GAGGT AGGGCC ACCT  GCATTTA-3’,  Reverse  (located  in  exon  1 1): 
5’-ATCCAATGCCTCGTAACAACCTGC-3’;  3)  for  exon  12  skipping  transcript:  Forward 
(located  at  junction  of  exon  1 1  and  exon  13):  5’-ATCTTAGTGGGCACAATAAAAG-3'; 
Reverse  (located  in  exon  13):  5’-CAGAAATTCTTGA  CCAGGTGCGGTA-3’.  The  relative 
intensity  of  the  PCR  products  was  measured  by  using  GeneTool  software  (Syngene). 

Testing  for  the  functional  significance  of  exclusion  of  BRCA2  exon  12  in  mouse 
embryonic  stem  (ES)  cells 

To  test  the  functional  significance  of  p.I2285V  and  BRCA2AE12  variants,  mutant  BAC 
DNAs  were  generated  in  the  human  BRCA2  BAC  RP1 1-777119  by  recombineering.  A 
mutant  BAC  DNA  containing  the  p.I2285V  variant  was  constructed  by  two  step  ‘hit  and 
fix’  method  as  described  previously  (Yang  and  Sharan,  2003).  The  BRCA2AE12  variant 
was  generated  in  the  BAC  by  using  galK  selection  and  counter-selection  method  as 
previously  described  (Warming  et  ah,  2005)  The  sequences  of  the  oligonucleotides  are 
available  upon  request.  Mutant  BAC  DNAs  were  then  introduced  in  PL2F7  mouse  ES  cells 
followed  by  deletion  of  the  endogenous  mouse  Brca2  gene  by  Cre-mediated  recombination 
(Kuznetsov  et  ah,  2008).  Recombinant  clones  lacking  endogenous  Brca2  were  selected  in 
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HAT  media  as  described  previously  (Kuznetsov  et  al,  2008).  We  tested  the  efficiency  of 
the  DNA  repair  function  by  challenging  the  recombinant  cells  with  DNA-damaging  agents: 
N-methyl-N’-nito-N-Nitrosoguanidine  (MNNG),  mitomycin  C  (MMC),  methylmethane 
sulfonate  (MMS),  cisplatin,  ultra-violet  light  and  y-irradiation.  Significance  of  the  survival 
difference  between  mutant  and  wild-type  cells  was  assessed  by  two-tailed  t-test  at  drug 
concentrations  corresponding  to  LD5o.  We  tested  the  efficiency  of  homologous 
recombination,  RAD5 1  foci  fonnation  and  karyotype  of  the  wild-type  and  the  recombinant 
cells  as  described  previously  (Kuznetsov  et  ah,  2008). 

For  the  RT-PCR,  we  synthesized  cDNA  using  RNA  isolated  from  mouse  ES  cells 
expressing  wild-type,  I2285V  and  Aexonl2  mutant  human  BRCA2  gene  using  Superscript 
III  reverse  transcriptase  (Invitrogen).  Region  flanking  exon  12  was  amplified  with  a 
forward  primer  B2exl  1FRT  (5’-CCAAGTCATGCCACACATTC-3’)  and  a  reverse  primer 
B2exl4RRT  (5’-ATTCTTGACCAGGTGCGGTA-3’)  from  exons  11  and  14,  respectively, 
and  the  PCR  conditions  included:  one  cycle  of  94°C  for  5  min,  35  cycles  of  94°C  for  30 
sec,  55°C  for  30  sec,  72°C  for  30  sec,  followed  by  incubation  at  72°C  for  5  min. 
Amplification  products  were  separated  in  a  1.6%  agarose  gel. 

Allelic  association  study 

As  controls,  the  study  enrolled  475  healthy  women  collected  through  the  New  York  Cancer 
Project  (NYCP).  The  NYCP  is  a  cohort  study  involving  consent  for  biospecimen  collection 
and  follow-up  of  8,000  healthy  volunteers  in  the  same  geographical  region  as  the  cases 


John  Wiley  &  Sons,  Inc. 


7 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Human  Mutation 


used  in  this  study  (Mitchell  et  ah,  2004).  For  the  cases,  we  genotyped  722  consecutive 
breast  cancer  cases  unselected  for  family  history  or  BRCA1/2  status  seen  at  Memorial 
Sloan-Kettering  Cancer  Center  and  collected  for  the  study  under  IRB  approved  protocols. 
All  individuals  in  the  study  were  self  identified  as  Ashkenazi  Jewish.  For  genotyping  of 
BRCA2  I2285V  we  used  TaqMan  allelic  discrimination  assay  under  standard 
manufacture’s  protocol  (Applied  Biosystems).  The  alleles  were  called  using  SDS  2.1 
software  (Applied  Biosystems). 

Results 

Conservation  and  co-occurrence 

When  we  examined  the  amino  acid  sequence  at  position  p.2285  of  human  BRCA2  with 
Align-GVGD  (http://agvgd.iarc.fr/alignnients/)  (Tavtigian  et  ah,  2006;  Mathe  et  ah,  2006), 
we  found  it  to  be  relatively  conserved  across  vertebrate  BRCA2  orthologs.  This  residue  is 
located  within  a  conserved  motif  encoded  by  exon  12.  From  humans  through  to  pufferfish 
(tetraodon),  this  position  is  either  lie  or  Leu.  Val  at  this  position  can  be  considered  to  be  a 
relatively  conservative  substitution  because  the  three  amino  acids  are  structurally  similar, 
each  with  a  nonpolar  side  chain. 

We  confinned  that  c.6853A>G  co-occurs  with  a  deleterious  variant  by  testing  85 
Ashkenazi  Jewish  individuals  from  Montreal  and  Toronto  who  had  previously  been  found 
to  carry  the  61 74delT  (c.5946delT)  variant  in  BRCA2.  We  identified  one  living  person  who 
carries  both  variants  (supplementary  figure  1).  We  re-contacted  this  person  and  established 
a  lymphoblastoid  cell  line  (LCL),  hereon  referred  to  as  DBV  (for  double  variant).  Being 
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unable  to  study  other  family  members,  we  could  not  establish  phase  genetically.  Instead,  by 
using  RNA  from  this  cell  line,  we  showed  that  exon  12  inclusion  and  exclusion  transcripts 
were  produced  from  both  alleles,  but  the  degree  of  exclusion  was  much  greater  with  for  the 
I2285V-associated  allele  than  for  the  c.6174delT- associated  allele.  We  were  able  to 
demonstrate  that  these  variants  occur  in  trans  (supplementary  Figure  2).  This  supports  the 
previous  segregation  and  family  history  data  showing  that  p.I2285V  is  very  unlikely  to  be  a 
highly  penetrant  allele  (as  discussed  above).  Nevertheless,  we  were  left  with  the  question  of 
the  splicing  out  of  exon  12  -  is  it  frequently  associated  with  BRCA2:p.I2285V,  and  if  so, 
what  is  the  consequence  of  this? 

Effect  of 1228 5 Von  the  BRCA2  mRNA 

The  clinical  protein  truncation  test  (PTT)  of  BRCA2  on  one  individual  carrying  this  variant 
indicated  that  exon  12  skipping  occurs  to  a  greater  degree  in  this  individual  than  the  level 
in  the  control  samples  (data  not  shown).  We  searched  for  splicing  regulatory  elements 
located  in  exon  12  by  using  the  following  on-line  prediction  tools:  ESEFinder: 
http://rulai.cshl.edu/tools/ESE:  RESCUE-ESE:  http://genes.mit.edu/burgelab/rescue-ese/ 
ESRsearch:http://ast.bioinfo.tau.ac.il/ESR.htm  and 

PESX:  http://cubweb.biology.columbia.edu/pesx/.  Using  RESCUE-ESE,  we  identified 
putative  ESEs  at  the  region  where  the  BRCA2:  c.6853A>G  variant  is  located.  Notably,  the 
variant  “G”  allele  disrupts  three  of  overlapping  putative  ESEs  (Figure  1A).  By  studying 
exon  12  retention  in  4  individuals  who  carry  the  p.I2285V  variant  and  two  controls,  we 
were  able  to  show  quantitatively  that  this  variant  significantly  increases  the  exclusion  of 
exon  12  of  BRCA2  (Figure  IB).  The  test  we  used  here  was  to  see  if  the  variant  was  retained 
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in  exon  12  by  capillary  sequencing  of  exon  12  inclusion  transcripts,  as  a  variant  promoting 
exon  12  skipping  would  result  in  a  diminished  signal  of  itself  in  the  retained  exon  12.  As 
shown  in  Figure  IB,  among  3  heterozygotes  (I2285V-a,  I2285V-b  and  I2285V-c),  the 
c.6853A>G  variant  (resulting  in  p. 12285 V)  is  barely  detectable  in  retained  exon  12, 
indicating  the  “G”  variant  strongly  promotes  exon  12  skipping.  This  causes  extremely 
skewed  A:G  allele  ratio  in  retained  exon  12,  to  a  degree  that  the  G  allele  (C  as  sequenced 
with  a  reverse  primer)  is  below  the  detection  power  of  capillary  sequencing.  The 
interpretation  of  the  results  is  somewhat  more  complicated  for  the  lymphoblastoid  cell  line 
DBV  (Figure  IB,  fifth  sample),  as  this  cell  line  is  derived  from  individual  compound 
heterozygote  for  c.6853A>G  and  c.6174delT  (discussed  above  and  shown  in  supplementary 
figure  2).  It  is  possible  that  the  ratio  of  A:G  has  been  reverted  by  the  presence  of  a  degree 
of  nonsense-mediated  decay  of  the  A  allele  intact  transcript,  resulting  from  the  upstream 
c.6174de/T  mutation  on  this  allele,  which  introduces  a  premature  stop  codon.  Despite  the 
fact  that  the  upstream  c.  61 74delT  may  equally  causes  some  decay  of  the  BRCA2AE12 
transcript  produced  from  this  allele  in  individual  DBV  (in  which  c.6853A>G  and 
c.6174delT  are  present  in  trans),  a  greater  degree  of  exon  12  skipping  was  observed  in  this 
cell  line  (as  the  sum  of  the  skipping  products  from  both  alleles)  than  in  the  wild-type 
counterpart  lymphoblastoid  cell  line  WLT  (Figure  1C),  suggesting  that  c.6853A>G 
enhances  exon  12  skipping.  Therefore,  the  biologically  relevant  outcome  of  this  single 
base  substitution  is  likely  to  be  the  disruption  of  a  putative  ESE,  resulting  in  increased 
exclusion  of  exon  12.  The  creation  of  a  novel  amino  acid  could  also  have  some  biological 
significance,  but  this  seems  less  likely  as  the  change  of  I  to  V  is  only  just  outside  the 
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tolerated  range.  Our  next  goal  was  to  establish  if  this  splicing  out  of  exon  12  has  any 
functional  consequences. 

Functional  test  of  c.6853A>G  variant  and  complete  deletion  of  exon  12 
To  test  the  functional  significance  of  deletion  of  exon  12,  mouse  ES  cell  based  functional 
assay  was  used.  The  B R CA 2 AE 1 2  variant  was  generated  in  the  BAG  DNA  coding  for 
human  BRCA2  and  introduced  into  mouse  ES  cells  containing  a  single  functional  copy  of 
Brca2  that  can  be  deleted  by  Cre-mediated  recombination.  Mouse  ES  cell  are  not  viable 
without  a  functional  copy  of  Brea 2.  After  Cre-mediated  deletion  of  endogenous  mouse 
Brca2,  we  tested  the  ability  of  human  transgene  to  rescue  the  lethality  of  Brca2- null  ES 
cells.  We  tested  6  independent  ES  cell  clones  expressing  the  BRCA2 AE12  transgene  and  5 
clones  expressing  c.6853A>G  variant  (resulting  in  p.I2285V)  for  viability  and  found  no 
difference  compared  to  the  wild-type  (data  not  shown).  RT-PCR  analysis  of  mouse  ES  cells 
expressing  BRCA2 AE12  variant  confirmed  the  deletion  of  96  bp  encoded  by  exon  12 
(Figure  2,  lanes  5  and  6).  These  cells  do  not  express  the  wild-type  transcript.  Similarly,  ES 
cells  expressing  the  p.I2285V  variant  resulted  predominantly  in  transcripts  lacking  exon  12 
(Figure  2,  lanes  3  and  4,  lower  band).  However,  this  variant  also  resulted  in  reduced 
amount  of  transcripts  that  retained  the  exon  12  (Figure  2,  lanes  3  and  4,  upper  band).  (Note 
the  presence  of  some  B RCA 2 AE  1 2  product,  mirroring  the  situation  in  vivo  in  humans). 
Taken  together,  these  results  show  that  the  deletion  of  exon  12  of  BRCA2  does  not  affect 
the  viability  of  mouse  ES  cells.  We  then  tested  the  Brca2- null  ES  cells  expressing  the 
BRCA2 AE12  variant  for  other  known  functions  of  BRCA2.  In  ah  experiments,  we 
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compared  the  results  for  the  B RCA 2 AE 1 2  variant  with  those  obtained  for  clones  that 
expressed  p.I2285V  and  wild-type  BRCA2. 

Two  clones  were  then  tested  extensively  for  sensitivity  to  DNA-damaging  agents  and  we 
found  no  significant  difference  for  either  the  BRCA2AE12  allele  or  the  p.I2285V  allele 
compared  to  the  control  cells  for  any  of  the  six  agents  used  (Figure  3,  panels  A-F).  As 
BRCA2  has  a  role  in  the  recruitment  of  RAD5 1  to  the  DNA  repair  sites  (Yuan  et  ah,  1999), 
we  tested  the  effect  of  the  deletion  of  exonl2  of  BRCA2  and  of  the  p.I2285V  allele  on 
radiation  induced  RAD5 1  foci  fonnation.  Both  the  BRCA2AE12  and  the  p. 12285  V  allele- 
containing  cells  formed  similar  a  number  of  foci  as  the  cells  expressing  wild-type  BRCA2 
(Figure  4).  The  known  role  of  BRCA2  in  homologous  recombination  and  genomic  stability 
(Khanna  and  Jackson,  2001)  prompted  us  to  test  these  cells  for  their  efficiency  of 
homologous  recombination  and  genomic  instabilities.  For  the  former,  we  examined  the 
efficiency  of  homologous  recombination  using  a  gene-targeting  vector  as  described 
previously  (Kuznetsov  et  ah,  2008).  This  targeting  vector  targets  Gt(ROSA)26Sor 
(abbreviated  here  as  Rosa26)  locus.  We  observed  similar  targeting  efficiency  in  the  ES  cells 
expressing  BRCA2AE12,  p.I2285V  or  wild-type  BRCA2  (Figure  5  A).  The  karyotyping  of 
these  cells  also  did  not  reveal  any  significant  increase  of  chromosomal  abnonnalities  due  to 
deletion  of  exonl2  (Figure  5B  and  C).  Based  on  these  results  we  conclude  that  the  deletion 
of  exonl2  of  BRCA2  does  not  significantly  affect  any  of  its  known  functions.  We  have 
previously  demonstrated  the  utility  of  the  ES  cell  based  assay  to  detect  subtle  functional 
defects  in  variants  that  are  predicted  to  be  of  low  risk  (Kuznetsov  et  ah,  2008).  However, 
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the  possibility  of  clinically  relevant  effects  of  this  allele  on  breast  cancer  risk  that  were  not 
due  to  the  known  functions  of  BRCA2  that  we  assayed  here  still  remained,  so  we  probed 
the  function  of  p.  12285V  from  a  clinical  perspective. 

Allelic  association  study 

The  loss  of  exon  12  could  have  an  effect  on  breast  cancer  risk  that  might  be  insufficiently 
large  to  be  evident  from  pedigree  analysis,  but  might  still  be  clinically  relevant.  To  address 
this  question,  we  assessed  the  frequency  of  BRCA2:  c.6853A>G  in  on-going  series  of 
Ashkenazi  Jewish  breast  cancer  cases  and  ethnically  matched  controls  (total  n  =  1 197).  We 
identified  10  carriers  in  475  controls,  and  15  carriers  in  722  cases,  resulting  in  an  odds  ratio 
for  the  association  between  this  allele  and  breast  cancer  of  0.99,  95%  confidence  interval 
0.44-2.21,  P  =  0.97.  When  we  adjusted  for  age,  the  odds  ratio  became  0.97  (95%  Cl:  0.41- 
2. 12),  P  =  0.88.  Setting  the  alpha  at  5%,  a  study  of  this  size  has  an  80%  power  to  rule  out 
an  allele  frequency  of  0.053  or  greater  in  the  cases,  when  the  allele  frequency  in  controls  is 
0.021. 

Discussion 

Here,  we  describe  the  analysis  of  BRCA2:  c.6853A>G  (BIC  nomenclature,  7081A>G, 
p.  12285V),  a  previously  reported  VUS.  We  first  confirmed  that  BRCA2:  c.6853A>G  co¬ 
occurs  in  trans  with  other  deleterious  alleles  in  BRCA2.  We  also  showed  that  this  variant 
promotes  increased  exclusion  of  exon  12  of  BRCA2,  but  is  not  associated  with  a  clinically 
relevant  increase  in  breast  cancer  risk.  Importantly,  we  have  carried  out  a  comprehensive 
analysis  of  the  functional  consequences  of  loss  of  exon  12  of  BRCA2,  using  a  murine  ES 
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cell  system,  wherein  functional  BRCA2  is  required  to  rescue  Brca2-null  murine  ES  cells 
(Kuznetsov  et  al.,  2008). 

Inappropriate  splicing  out  of  functionally  critical  exons  is  associated  with  disease.  Most  of 
the  mutations  that  affect  splicing  are  those  that  disrupt  splice  sites  (Hastings  and  Krainer, 
2001).  More  recently,  however,  exonic  mutations  that  affect  splicing  have  been  identified 
(Cartegni  et  ah,  2002).  Our  results  show  that  BRCA2 :  c.6853A>G  is  one  such  variant,  and 
that  a  degree  of  exclusion  of  exon  12  is  seen  in  alleles  without  this  variant  (Figures  1  and 
2).  It  appears  that  ESEs  may  be  present  in  exon  12  of  BRCA2  (Figure  1A),  and  it  is  notable 
that  12285,  and  the  surrounding  amino  acids  (KRRGEPS  I  KRNLLNEFDR)  are  well 
conserved  through  vertebrate  evolution.  Thus,  the  question  of  the  significance  of  loss  of 
exon  12  is  an  interesting  one.  As  noted,  loss  of  exon  12  results  in  a  BRCA2  protein  that  is 
shorter  than  the  wild-type  form  only  by  32  amino  acids,  and  this  part  of  BRCA2  does  not 
contain  any  known  important  functional  motifs.  Nevertheless,  BRCA2AE12  has  been 
reported  to  be  more  commonly  found  in  breast  cancers  than  in  normal  tissues  from  the 
same  individuals  (Bieche  and  Lidereau,  1999).  In  this  study  of  breast  cancer  and  matched 
normal  tissue  from  12  women,  the  proportion  of  BRCA2AE12  transcripts  was  higher  in  the 
tumor  in  four  pairs,  and  in  no  pairs  was  the  BRCA2AE12  higher  in  the  normal  tissue  than 
in  the  tumor  tissue.  Similarly,  in  a  study  of  7  ovarian  cancer  cell  lines,  5  were  found  to  have 
a  two-  to  five-fold  increase  in  BRCA2AE12  compared  with  mean  values  found  in  cultures 
of  normal  human  ovarian  surface  epithelium.  The  same  over-expression  in  tumor  cell  lines, 
rather  than  nonnal  epithelium,  was  seen  in  one  of  four  breast  cancer  cell  lines  and  one  of 
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four  prostate  cancer  cell  lines.  In  the  case  of  DU-145,  the  wild-type  isoform  of  BRCA2  was 
barely  detectable,  whereas  the  truncated  isofonn,  BRCA2AE12,  was  strikingly  over¬ 
expressed  (Rauh-Adehnann  et  ah,  2000).  Despite  these  findings,  the  functional  analyses 
presented  here  indicate  that  BRCA2AE12  can  assume  the  biological  role  of  full-length 
BRCA2.  Very  likely,  the  increased  level  of  BRCA2AE12  in  cancer  cells  is  an 
epiphenomenon  rather  than  a  causative  mechanism. 

We  do  not  think  BRCA2AE12  represents  a  clinically  important  splice  variant  of  BRCA2. 
Nevertheless,  it  is  important  to  investigate  all  potentially  deleterious  intra-exonic  VUS 
using  several  approaches.  For  example,  the  BRCA2  variant  c.8393C>G  (p.T2722R)  was 
initially  found  to  result  in  an  out-of- frame  fusion  of  exons  17  and  19  of  BRCA2  (Fackenthal 
et  ah,  2002).  Later,  it  was  found  that  the  skipping  of  exon  18  is  not  complete  (erratum,  Am 
J  Hum  Genet,  December  2003)  and  thus  the  significance  of  the  original  finding  was 
questioned.  Subsequently,  both  genetic  (Easton  et  ah,  2007)  and  functional  (Kuznetsov  et 
ah,  2008)  analyses  have  confirmed  that  this  variant  is  deleterious. 

Although  the  allelic  association  study  we  have  carried  out  has  not  ruled  out  a  small  effect 
on  risk,  14,199  Ashkenazi  Jewish  cases  and  an  equal  number  of  controls  would  be  required 
to  rule  out  an  odds  ratio  of  1.25  (with  the  same  1-P  and  a  as  above).  Therefore,  it  remains 
possible  that  the  splicing  out  of  exon  12  does  result  in  some  subtle  perturbations  in  risk  for 
breast  cancer,  but  at  the  current  time,  we  cannot  detect  such  an  effect,  if  is  exists  at  all. 


John  Wiley  &  Sons,  Inc. 


15 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Human  Mutation 


The  Breast  Cancer  Information  Core  (BIC)  (http://research.nhgri.nih.gov/bic/)  provides 
easily  accessible  infonnation  on  numerous  variants,  but  is  dependent  upon  timely  entry  of 
new  information  about  all  variants.  Based  on  a  new  classification  system  for  VUS  that  was 
determined  by  an  IARC  working  group  and  recently  published  (Plon  et  ah,  2008),  we 
suggest  that  this  variant  should  now  referred  to  a  class  1  variant  (probability  of  being 
pathogenic  of  less  than  0.1%).  This  conclusion  is  amply  supported  by  other  recently- 
published  model-based  evidence  (Spearman  et  ah,  2008). 

Here,  we  show  for  the  first  time,  that  loss  an  exon  of  a  tumor  suppressor  gene  has  no 
measurable  consequences  across  a  range  of  activities.  The  technique  described  here  could 
be  used  to  comprehensively  assess  the  separate  functionality  of  each  exon  of  BRCA2,  and 
the  findings  raise  the  question  of  a  possible  class  of  exons  that  are  dispensable  for  BRCA2 
tumor  suppressing  functions.  Indeed,  earlier  biochemical  studies  have  revealed  that 
heterologous  fusion  proteins  consisting  of  BRCA2  domains  interacting  with  Rad5 1  and 
Dssl  are  capable  of  restoring  DNA  repair  in  BRCA2  deficient  cells.  The  BRCA2  domains 
mediating  association  with  Rad5 1  and  Dssl  are  the  BRC  repeats  (central)  and  OB 
(oligosaccharide-binding,  C-terminal)  domains,  respectively  (Saeki  et  ah,  2006).  Notably, 
the  peptide  encoded  by  exon  12  is  not  involved  in  either  of  these  domains.  Thus  the  murine 
ES  system  used  here  could  be  used  to  investigate  which  functions,  if  any,  are  impaired  in 
cells  that  lack  various  combinations  of  BRCA2  exons.  It  is  interesting  to  note  that  exons  10 
and  12  of  BRCA2  do  not  encode  functionally  critical  peptides  and  are  in-frame  (i.e.  the 
reading  frame  will  not  be  altered  when  these  exons  are  not  represented  in  the  transcript). 
Therefore,  missense  variants  located  in  these  two  exons  that  alter  an  amino  acid  and/or 
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cause  skipping  of  the  entire  exon  are  likely  neutral.  This  statement  should  be  tempered  by 
the  possibility  that  variants  could  exist  that  may  have  unpredictable  effect  on  the  stability  or 
structural  integrity  of  BRCA2  protein.  However,  such  variants  are  likely  to  be  rare.  In 
support  of  this  hypothesis,  on  reviewing  the  published  literature  (Easton  et  ah,  2007; 
Spearman  et  ah,  2008;  Gomez  Garcia  et  ah,  2009)  and  the  BIC 
http://research.nhgri.nih.gov/bic/,  we  were  unable  to  identify  any  known  pathogenic 
misssense  variants  in  exons  10  or  12  of  BRCA2.  While  not  conclusive,  this  observation 
supports  our  contention  that  exon  12  (and  possibly  exon  10)  is  functionally  dispensable. 
Using  the  murine  ES  cell  system  reported  here,  it  will  be  interesting  to  see,  if  such 
redundancy  is  seen  in  other  multi-exon  tumor  suppressor  genes. 
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Figure  Legends 

Figure  1:  Splicing  at  BRCA2:  c.6853A>G  in  humans 

A)  Potential  exonic  splicing  enhancer  (ESEs)  disrupted  by  BRCA2:  c.6853A>G 
(p.I2285V) 

RESCUE-ESE  (http://genes.mit.edu/burgelab/-rescue-ese/)  identified  putative  ESEs  at  the 
region  where  the  BRCA2:  c.6853A>G  variant  is  located.  Notably,  the  variant  “G”  allele 
disrupts  three  of  overlapping  putative  ESEs. 

B)  Splicing  of  BRCA2  exon  12  in  humans  carrying  the  c.6853A>G  allele 

The  amplicon  of  exon  12-13  was  sequenced  using  a  reverse  primer  located  at  exon  13.  The 
position  of  c.6853  in  the  sequencing  chromatogram  is  indicated  by  a  dashed  box.  The  wild- 
type  A  allele  (shown  here  as  T  on  the  reverse  strand)  is  the  major  signal  in  the  three 
c.6853A>G  carriers  (I2285V-a,  I2285V-b  and  I2285V-c)  and  a  non-carrier  control 
(control).  The  cell  line  DBV  carries  double  variants  of  6174delT  and  c.6853A>G.  The 
result  shown  for  DBV  suggests  that  there  is  no  technical  allele  bias  in  the  PCR.  WLT  is  a 
wild-type  control  line  carrying  neither  variant. 

C)  Relative  quantification  of  exon  12  skipping  in  WLT  and  DBV  LCLs 

Three  amplicons  were  produced:  El 2- 13  with  exon  12  inclusion,  AE12  with  exon  12 
skipping  and  El  1  as  a  control  of  endogenous  expression  level.  The  primers  are  listed  in  the 
methods  section.  C:  Semiquantitative  comparison  of  exon  12  skipping  in  WTL  and  DBV: 
normalization  was  done  by  considering  Ell  expression  as  1.00.  The  relative  expression 
level  of  E12  intact  vs.  AE12  is  shown  as  a  percentage  bar.  Blue  bars:  exon  12  skipping 
transcripts,  black  bars:  exon  12  intact  transcripts.  The  error  bar  indicates  the  variance  seen 
following  four  measurements. 
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Figure  2:  RT-PCR  demonstrating  96-bp  deletion  in p.  I2285V and BRCA2AE12  mutants. 

Ethidium  Bromide  stained  gel  showing  RT-PCR  results  of  two  mouse  embryonic  stem  cell 
clones  generated  independently  for  each  variant.  WT :  wild-type,  negative  control:  water 
control,  FL:  fragment  of  a  full  length  product,  AE12:  transcript  with  deletion  of  exon  12. 
Numbers  below  the  gel  picture  are  lane  numbers.  Marker  is  1  kb  plus  ladder  (Invitrogen). 

Figure  3:  Sensitivity  of  ES  cells  expressing  mutant  and  wild-type  BRCA2  to  different 
DNA-damaging  agents. 

Survival  of  ES  cells  expressing  p.I2285V,  BRCA2AE12  mutants  compared  to  ES  cells 
expressing  wild  type  BRCA2  after  exposure  to  different  DNA-damaging  agents:  Methyl- 
N’-nitro-N-Nitrosoguanidine  (MNNG,  A),  mitomycin  C  (MMC,  B),  cisplatin  (C),  methyl- 
methanesulfonate  (MMS,  D),  ionizing  radiation  (IR,  E)  and  ultraviolet  light  (UV,  F).  Drug 
sensitivity  was  expressed  as  a  percentage  of  surviving  cells  compared  to  untreated  cultures. 
Two  independently  generated  clones  for  each  variant  were  tested  and  their  average  values 
are  shown  here.  Error  bars  indicate  standard  deviations.  P  values  were  equal  0.19693 
(I2285V);  0.16019  (BRCA2AE12)  for  5  pM  MNNG,  0.09178  (I2285V);  0.08212 
(BRCA2AE12)  for  10  ng/ml  MMC,  0.5494  (I2285V);  0.41 10  (BRCA2AE12)  for  0.4  pM 
cisplatin,  0.1 1846  (I2285V);  0.74637  (BRCA2AE12)  for  15  pg/ml  MMS,  0.1 131  (I2285V); 
0.14066  (BRCA2AE12)  for  200  Rad  IR  and  0.3964  (I2285V);  0.4238  (BRCA2AE12)  for 
10J/m2  UV. 
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Figure  4:  Radiation-induced  RAD51  foci  formation  in  mouse  ES  cells  expressing 
BRCA2  transgenes. 

A)  RAD51  foci  formation  before  (-IR)  and  after  (+IR)  ionizing  radiation.  RAD5 1  foci  are 
shown  in  green,  y-H2AX  foci  marking  sites  of  DNA  damage  are  shown  in  red,  nuclei  are 
stained  with  DAPI  (blue).  The  right-most  panel  shows  the  merged  picture.  Two 
independent  clones  showed  similar  results  and  only  one  of  them  is  shown  in  each  case. 

B)  Quantification  of  RAD51  and  y-H2AX  foci  before  (-IR)  and  after  (+IR)  ionizing 
radiation.  30  nuclei  were  counted  in  each  case  and  their  mean  values  are  shown.  Error  bars 
indicate  means  ±  s.d. 


Figure  5:  Effect  of  BRCA2  variants  on  homologous  recombination  and  genomic 
integrity  in  mouse  ES  cells. 

a)  Homologous  recombination  efficiency  shown  by  gene  targeting  at  Rosa26  locus. 
Numbers  above  the  bars  indicate  the  actual  number  of  colonies  undergoing  homologous 
recombination  compared  to  total  number  of  colonies  tested. 

b)  Karyotype  analysis  of  p.I2285V,  BRCA2AE12  mutants  compared  to  WT  control  cells. 
Representative  metaphases  of  the  mutants  and  WT  control  cells  are  shown.  Ten  metaphases 
were  counted  in  each  case  and  scored  for  different  chromosomal  abnonnalities.  c)  The 
identified  chromosomal  aberrations  are  tabulated. 
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Figure  1:  Splicing  at  BRCA2:  c.6853A>G  in  humans 
A)  Potential  exonic  splicing  enhancer  (ESEs)  disrupted  by  BRCA2:  c.6853A>G  (p.I2285V) 
B)  Splicing  of  BRCA2  exon  12  in  humans  carrying  the  c.6853A>G  allele 
C)  Relative  quantification  of  exon  12  skipping  in  WLT  and  DBV  LCLs 
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Figure  2:  RT-PCR  demonstrating  96-bp  deletion  in  p.  I2285V  and  BRCA2AE12  mutants. 
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Figure  3:  Sensitivity  of  ES  cells  expressing  mutant  and  wild-type  BRCA2  to  different  DNA-damaging 

agents. 
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Figure  4:  Radiation-induced  RAD51  foci  formation  in  mouse  ES  cells  expressing  BRCA2  transgenes. 
A)  RAD51  foci  formation  before  (-IR)  and  after  (+IR)  ionizing  radiation. 

B)  Quantification  of  RAD51  and  y-H2AX  foci  before  (-IR)  and  after  (+IR)  ionizing 
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Figure  5:  Effect  of  BRCA2  variants  on  homologous  recombination  and  genomic  integrity  in  mouse  ES 

cells. 

a)  Homologous  recombination  efficiency 

b)  Karyotype  analysis  of  p.I2285V,  BRCA2AE12  mutants  compared  to  WT  control  cells. 
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Supplementary  Figure  legends 
Supplementary  Figure  SI 

Pedigrees  of  four  families  carrying  BRCA2:  p.  I2285V 

Figure  la 

Individual  V :2  presented  to  the SM  BD -Jewish  General  Hospital,  M  ontreal,  Quebec,  requesting 
genetic  testing  for  a  BRCA2  variant  that  he  said  has  been  identified  in  his  family.  He  brought 
along  genetic  testing  results  which  allowed  us  to  identify  the  variant  as  12285V.  The  family  is  of 
Ashkenazi  Jewish  origins  and  was  originally  tested  in  another  Canadian  province.  Individuals 
V  :3  and  IV  :2  had  already  been  found  to  carry  this  variant.  We  declined  to  offer  clinical  testing  to 
V:1  but  asked  IV:  2  to  come  in  for  more  detailed  studies.  Having  confirmed  the  presence  of  the 
variant,  samples  from  this  individual  were  used  to  further  characterize  the  variant. 

Figure  lb 

Individual  111:1  presented  for  genetic  testing  at  the  M  cGill  University  Health  Centre,  M  ontreal, 
Quebec,  on  account  of  her  diagnosis  of  serous  papillary  ovarian  cancer  diagnosed  at  age  44, 
together  with  her  Ashkenazi  Jewish  ancestry.  A  three  generation  pedigree  was  obtained  and 
pathology  records  confirmed  the  diagnosis  in  111:1.  Pathology  records  were  not  available  for  the 
remaining  relatives.  Genetic  testing  for  the  three  Ashkenazi  Jewish  founder  mutations  was 
offered  to  individual  111:1  and  revealed  the  presence  of  a  B  RCA 2:6174delT  mutation.  Following 
identification  of  a  BRCA2  mutation  in  the  family,  individual  111:2  presented  for  genetic  testing. 
Subsequent  testing  for  the  three  founder  mutations  also  identified  individual  111:2  as  a  carrier  of 
the  B  RCA  2  6174delT  mutation.  At  the  time  of  testing,  111:2  also  provided  a  blood  sample  for 
research  purposes.  In  the  course  of  our  efforts  to  identify  a  compound  heterozygote,  single-stand 
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conformation  analysis  revealed  the  presence  of  the  B  RCA 2  12285V  variant.  Further  testing  of 
the  lymphocyte  DN  A  from  individual  111:1  also  revealed  the  existence  of  the  B  RCA 2  12285V 
variant.  No  other  family  members  were  available  for  testing. 

Figures  lc  and  Id 

Both  families  all  have  Ashkenazi  Jewish  ancestry  and  all  had  full  sequencing  of  B  RCA  1/2  by 
M  yriad  Genetics  Laboratories  (Salt  Lake  City,  UT).  In  Figure  lc,  individual  lll.l  was  seen  for  an 
initial  visitand  tested  in  2004,  when  the  BRCA2: 12285V  variant  was  identified.  In  Figure  Id, 
lll.l  was  tested  elsewhere  in  2007  and  the  BRCA2:  12285V  variantwas  identified.  Later  that 
year,  the  proband  came  to  M  emorial  Sloan  Kettering  Cancer  Center  for  an  initial  visit  to  discuss 
this  result. 

NB  in  all  figures  the  proband  is  indicated  with  an  arrow. 

Supplementary  Figure  S2 

A  double  variant  lymphoblastoid  cell  line  (LCL)  referred  to  as  DBV  was  established  from  an 
individual  carrying  both  BRCA2:  NM_000059.3:c.5946delT  (p.Serl982A rgfsX 20,  thefounder 
mutation  known  as  6174delT,  according  to  BIC)  and  BRCA2:  NM_000059.3:  c.6853A>G 
(p. 12285V,  7081A>G  according  to  BIC).  The  top  two  rows  of  sequencing  electrograms  showed 
the  presence  of  the  double  variants  in  the  genomic  DNA  of  this  cell  line.  For  genomic  DNA 
sequencing  ,  primer  pairs  of  forward  5'-AGGTTGTTACGAGGCATTGG-3'  with  reverse  5'- 
TCCAGAGAA  AG  CAGATG  AATTT-31;  and  forward  51- 

TCTGTGGTATCTGGTAGCATCTG-31  with  reverse  5'-GCACA  GTGGCTCATGTCTGT-31 , 
were  used  to  amplify  the  DNA  regions  of  c.5946delT  and  c.6853A>G,  respectively.  Forward 
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primers  were  used  as  sequencing  primers.  Allele-specific  PCR  with  cDNA  was  performed  with 
forward  primers  differentiating  thee. 5946T  wild-type  allele  and  c.5946delT  allele:  5'- 
GGGATTTTTAGCACAGCAAGG  -3'  (for  c.5946delT);  5'- 

TTGTGGGATTTTTAGCACAGCAAGT-3'  (for  the  wild  type  allele  at  C.5946T),  the  reverse 
primer  of  5'-GTGCGAAAGGGTACACAGGT-3'  was  used  for  amplifying  both  alleles. 

C.5946T  wild  type  allele-specific  PCR  was  performed  with  cDNA  from  a  DBV  LCL  and  a  wild- 
type  LCL,  both  exon  12  inclusion  (of  1085bp)  and  skipping  (of  989bp)  transcripts  were  detected. 
The  schema  for  the  allele-specific  PCR  reactions  is  shown  in  the  next  two  figures  (top, 
c.5946delT  allele;  bottom  c.6853A>G  allele). 

The  next  panel  shows  that  the  degree  of  exon  12  skipping  from  the  wild-type  allele  (C.5946T)  is 
higher  in  the  DBV  line  than  in  wild-type  individuals.  In  the  final  row,  cDNA  from  DBV  LCL 
was  used  for  allele-specific  PCR,  by  using  C.5946T  (normal)  and  c.5946delT  primer, 
respectively.  Sequencing  the  exon  12  inclusive  transcript  from  the  allele  of  c.5946delT  by  using 
the  reverse  primer  showed  "A"  (T  on  the  reverse  strand)  at  the  position  of  6853,  illustrating  that 
c.5946delT  andc.6853A>G  are  in  trans. 
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Figure  SI:  Pedigrees  of  four  independent  probands  carrying  the  BRCA2:  p. 12285V  allele 


Figure  la 
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Figure  1b 
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FigureS2:  BRCA2  founder  mutation  c.5946delT  and  variant  c.6853A>G  are  in-trans 
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